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Abstract 

The Faddeev Random Phase Approximation is a Green's function technique that makes use 
of Faddeev-equations to couple the motion of a single electron to the two-particle-one-hole and 
two-hole-one-particle excitations. This method goes beyond the frequently used third-order Alge- 
braic Diagrammatic Construction method: all diagrams involving the exchange of phonons in the 
particle-hole and particle-particle channel are retained, but the phonons are described at the level 
of the Random Phase Approximation. This paper presents the first results for diatomic molecules 
at equilibrium geometry. The behavior of the method in the dissociation limit is also investigated. 
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I. INTRODUCTION 



: first-principle calculations has taken a high 



-|5j. The Green's function formalism [6- 



The study of electronic systems by means o 
rise thanks to modern computer technology 

ij] is one of these first-principles methods that has been succesfully applied in quantum 
chemistry [9 12 ] . The correlations in a many-body system are described in terms of an 
electron self-energy that acts as an energy- dependent potential describing the motion of a 
single electron in the many-electron system. 

A particular third-order approximation scheme to the self-energy can be obtained using 
the Algebraic Diagrammatic Construction (ADC(3)) 13] method as developed by Schirmer 
and coworkers. This method has proven to be very successful in predicting one-electron 
properties in molecules [l4] as measured e.g. in electron momentum spectroscopy. Although 
the equations were derived in a purely algebraic manner, they can be shown to be equivalent 
to resumming all particle-hole (ph) and particle-particle (pp) interactions between two- 
particle-one-hole (2plh_) and two-hole-one-particle (2hlp) states up to the Tamm-Dancoff 
approximation (TDA) |7| level. This is diagrammatically equivalent to considering phonons 
(excitations in the ph and pp channel) at the TDA level, and then allowing the exchange of 
these phonons in all possible ways between the tree propagators describing the 2plh/2hlp 
states. 

The TDA allows no ground-state correlations in the construction of the phonons. An 
improvement in this respect is the Random Phase Approximation (RPA) I5j. Calculations 
for the electron gas show that these improvements lead to a correct prediction of the plasmon 
pole, whereas the TDA completely fails to describe the plasmon spectrum. It is therefore of 
interest to formulate an analogous theory to ADC(3) that resums the ph and pp interactions 
up to RPA level. 

Going beyond the TDA level has proven to be very difficult [if]], even though it is known 
that the RPA should be better to describe collective behavior, at least for nuclear sys- 
tems 15[. The Faddeev Random Phase Approximation (FRPA) [17| solves this problem by 
using the Faddeev technique to include RPA-phonons in the self-energy. This method has 



succesfully been applied to both nuclei 18 



191 ] and atoms [20(. It is the aim of the present 



paper to study the application of this technique to simple molecular systems. 

In the second section of this work we give a short overview of the working equations for 



the FRPA method. In section II III we present the numerical results for a set of diatomic 
molecules. A summary is provided in section IIVI 



II. THEORY 



A. Single-particle Green's Function 

The single-particle motion in an N-body system is described by the single-particle prop- 
agator 
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where T[...] represents the time-ordering operator, is the exact ground state and a a (t) 
and aj,(t) are the addition and removal operators in the Heisenberg representation for an 
electron in a single-particle state a. For practical calculations it is more convenient to use 
the Lehmann representation of the Green's function 
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where the ^ r ^ ±1 represent exact eigenstates of the Hamiltonian with energy E^ 1 . This 
transition to the energy domain transforms the Dyson equation from an integral equation 
into the algebraic relation 

Gaji (E) = G% (E) + £ ( E ) ( E ) G W ( E ) ■ 
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In this equation the exact Green's function G is expressed in terms of the non-interacting G^ 
and the irreducible self-energy Y<*(E). Approximation schemes for the single-particle Green's 
function boil down to finding an appropriate perturbation expansion for the irreducible self- 
energy. 

In our approach, we want to couple the single-particle states with 2plh and 2hlp states. 
According to Refs. 



21 



22j, the connection between the irreducible self-energy E* and the 
six-point response function R can be written as 

S a,/3 ( E ) = + T ^2 V ^M R ^u,e9a (E) V&fia (4) 



FIG. 1. The Feynman-diagram for the irreducible self-energy £* in equation ([5]) within the FRPA. 
The first diagram represents the HF-like static self-energy. 



where V is the anti-symmetrized two-particle interaction and S is the static self-energy 
as depicted in Figure [TJ We now replace the exact single-energy six-point response function 
R{E) by an approximate propagator that has indices that are restricted to the 2plh space 
(R p ) or 2hlp space (R 2hlp ), and that is exact up to third order: 

S a,/3 {E) = + 7 E E U aLV,\nR\p,v,<iOo (E) U e Q,p a . (5) 

The two-particle interaction V in Eq. fl4]) has been replaced by a second order expansion 

U a p t *fS — {\xp,\n + AU a p ; \n) Vxprf. (6) 

X,fi 

This AU is needed to guarantee full summation up to third order perturbation theory and 
was chosen to be the same as the vertex correction used in the ADC (3) 



B. pp/ph RPA interaction 



The two-particle propagator is denned by 
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E - (Z + ir) n E- en - 277 
where the A' pp , y pp and e pp are shorthand notations for the overlap amplitudes and energy 
differences in Eq. ([7]). A relevant approximation for this object is obtained by solving the 
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FIG. 2. The diagrammatical representation of the pp-RPA equation (a) and the ph-RPA equa- 
tion |(b)| where the single lines represent non-interacting and the double lines interacting propaga- 
tors. 
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as indicated diagrammatically in Figure 2(a) Equation ([TO]) defines the effective pp interac- 



tion T pp , which includes dynamical screening and will be used later as a building block for the 
2plh and 2hlp interaction. This simple form of the Bethe-Salpeter-like equation for the pp 
propagator in function of a screened interaction T pp is possible because the non-interacting 
pp propagator is diagonal in the HF basis. 

The same procedure can be followed for the particle-hole (ph) polarization propagator 



(see Figure 2(b)), denned as 

= E 





CtL<l a 


ml 






E — 


'E N - 

k m 




+ if] 



E 



E - W - E») 



E 



X ph X ph ^ yP h t yP h 



af3,m 7<5,m 



E - e*f + + i, 



b-E 



a/3,n^ 7<5,n 



z?7 



(11) 
(12) 



The corresponding Bethe-Salpeter-like equation in the RPA reads as 

(£) = n^S (£0 + £ n *S (*0 ^< 75 (E) 
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and defines the effective ph interaction T ph . 

The actual calculation of the amplitudes and poles of the pp propagator and ph polar- 



ization propagator can be done by solving the generalized eigenvalue problems 

a b \ f x pp+ ypp- \ ( 11 o \ / x pp+ ypp- 
B \ c ) \ ypp+ xpp- J \ o -li / \ y pp+ xpp- 



15| 




(15) 



where 



A a p n s = - S aS 5^) (e a + €p) + -V a p >y s a,P,j,5 > F (16) 

B a p n& =V a/3 ^ a, f3 > F; 7, 5 < F (17) 

Ca/3, 7 5 = (Sa-ySpS ~ baS^M) (^ + CjS) ~ ^Q/3, 7 5 Ot, /?, 7, 5 < F. (18) 

Here the e a represent Hartree-Fock single-particle energies with the Fermi level F separating 
the occupied and unoccupied HF levels. The equations for the ph polarization propagator 
are again very similar: 
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where 



D a p rf s =6crf6ps {e a - ep) + V aS ,^ a,j>F;/3,S<F (20) 

E a (3 n 8 =Vart,p8 <*, 7 > E ] ft , 5 < F. (21) 



C. Faddeev equations 

The diagrammatic content of R cannot be cast into the form of a Bethe-Salpeter equation 
without double counting of some classes of diagrams, in contrast to the more complicated 



4-times propagator (see Ref. llTJ). That is why the Faddeev technique [23] must be used to 
split this object into three parts. The analysis will be done for R 2 P lh (the derivation of R 2hlp 
is found to be completely analogous, but with an interchange of particle and hole lines). 
The decomposition of R 2plh into three Faddeev components R^ reads 

R2 a^Mu ( E ) = G aM,X^ ( E ) ~ G a/3y,fj,\v ( E ) + ^2 R *PlM» ^ ' ( 22 ) 

i=l,2,3 
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FIG. 3. Diagrammatic representation of equation (|25p . 



where (j(°) > is the part of the non-interacting 2plh propagator with positive energy 
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Together with its exchange counterpart, they form the free 2plh propagator 
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The relation between the different components can be derived from the diagrammatic 
content of Figure El The superscripts (i), (j) and (k) are cyclical permutations of 1, 2 and 
3 and correspond to the numbering of the fermion lines from left to right. In our notation 
lines 1 and 2 are the particles and line 3 is the hole. Each propagator RW ends with 
lines j and k interacting through the adequate RPA interaction vertex, while all possible 
prior propagation is included in Ry\ R^ and the non- interacting propagators. is the 
extension to 2plh space of T pp and T ph by adding a Kronecker delta for the third fermion 
line. The corresponding Bethe-Salpeter equations for the R^ 
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form a closed self-consistent system. 
The Lehmann representation 
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can be derived from the Lehman representation for the full R (see Ref. (17|). The sum of 
the R^f ree makes sure that the non-interacting poles appearing in the first term of Eq. ( 122]) 



are precisely cancelled. The spectroscopic amplitude can be recovered by summing over the 
three Faddeev components 

XnR*v m = N Xja • (27) 



i=l,2,3 



By multiplying equation (|2T)j) with (E — e™) and taking the limit for E — >■ e™, the problem 
is reduced to a non-linear eigenvalue problem for the spectroscopic amplitudes and the poles. 
The non-interacting poles do not coincide with the Faddeev-poles, so the Rf ree is guaranteed 
to disappear when taking the limit: 



(28) 



The explicit treatment of this equation for i = 3 (i.e. the pp channel) is given in appendix \A\ 
and is easily extended to the two other channels. When substituted in equation (|28|) . we 
arrive at 

X {i) = { u® _ 1 r (0t + #»#(i)t ) (^(i) + . (29) 



The vectors D«, T® and ff« are all diagonal in the freely propagating line and can 
be written in terms of the pp- and ph-amplitudes and energies. Their explicit form can be 
found in Ref. 17J. By introducing a vector containing these three components, 
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this non-linear equation in the Faddeev-energies and amplitudes can be written in the form 

(31) 
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where the matrix M, 
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takes care of the coupling between the different channels. After some matrix algebra, this 
can be converted into a linear non-hermitian eigenvalue problem 
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e Fd X = (11 - HH*M) 1 U [T ] M + DU~ l (ll - HH ] M)] X. (33) 
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The matrix dimension of the eigenvalue problem is three times the size of the 2plh-basis. 
Two thirds of the solutions are spurious and can be projected out, so the actual matrix 
dimension reduces to the size of a single 2plh basis. 



D. Handling spurious solutions 



The use of Faddeev-equations inherently introduces spurious solutions |24H26l| . The so- 
lutions for which the sum in Eq. (J27J) is zero, have no physical meaning and have to be 
discarded. At the same time the vectors themselves will have to be anti-symmetric under 
exchange of the two particle or hole lines. By projecting the Hamiltonian matrix ( 133]) onto 
the vector that has the right symmetry properties, and is non- vanishing when summed, the 
matrix dimension is reduced by a factor of 3. This vector space is spanned by the vector 
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San&px&iv The dimension of the matrix is now the same as in the 



where (ll e 

standard ADC(3) matrix problem 13j. It can be verified that by using Tamm-Dancoff 
(TDA) interactions and after performing this projection, one regains the ADC(3) equations 
(see Appendix IBl). 



E. Single particle propagator and ground-state properties 



The calculation of the FRPA single-particle propagator is now done by diagonalization 
of the symmetric matrix 



(35) 



where the e Fd matrices are diagonal and contain the 2plh and 2hlp Faddeev energies. The 
tilde indicates that the coupling matrix elements are written in the basis that diagonalizes 





p/h 


2plh 


2hlp 


p/h 


( e 


U 


V \ 


2plh 




e Fd 





2hlp 










Ua t m ^ ^ U QtvX^Xx^viYi. (*^) 



the Faddeev matrices: 



Note that in standard ADC(3) it is possible to write the equivalent of matrix fl35|) using (1B6|) 
and ([6]) as sub-blocks without a separate diagonalization in 2plh and 2hlp space. This is not 
the case in the FRPA formalism as due to the nonhermiticity of the right hand side of Eq. [33 
Thus, one should first diagonalize the 2plh and 2hlp sub-blocks (that is, solve the Faddeev 
equations) and then write the matrix (135]) in the new basis obtained. Performing the double 
diagonalization procedure therefore involves a doubling of the computer time with respect 
to the usual ADC(3) approach. In practical calculations, however, this is not the case since 
the dimension of matrix ( |35i) can be reduced drastically by employing Arnoldi techniques in 
the 2plh and 2hlp diagonalizations. This approach has been applied previously 19[ and it 
was found that a limited number of Arnoldi vectors guarantee correct converged results for 
total energies and ionizations potentials. In this paper, however, we did not resort to the 
Arnoldi algorithm and all results are obtained with full diagonalizations. 

The diagonalization of ( 1351) results in energies u n and residues f a , n (see Eq. (T5])), defining 
a new single-particle Green's function. By summing over the the solutions below the Fermi 
level, the density matrix 

n a ,p = /«.n/j9,n ( 37 ) 

n<F 

and the corresponding ground-state energy 

< = i I > J "nf a , n f* I 
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can be obtained. 

In principle full self-consistency could be achieved by iteratively recalculating the phonons 
on the basis of the new propagator and applying the Faddeev procedure. This is however 
computationally too demanding. We do improve the self-consistency of the solution by 
updating the Hartree-Fock-like static self-energy diagram. Instead of the diagonal matrix of 
single-particle energies, the Hartree-Fock self-energy calculated with the new density matrix 
n a fi has to be included in the diagonalization. Note that, both in FRPA and ADC(3), this 
partially self-consistent treatment is needed to include all static self-energy diagrams up to 
third order. 
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III. RESULTS AND DISCUSSION 



The accuracy of the FRPA method is evaluated by comparing to the ADC(3) method, 
and to Coupled Cluster calculations with Single, Double and Perturbative Triple [CCSD(T)] 
excitations. The latter method should be of a comparable level of theory as both the ADC(3) 



and FRPA. Where possible, the comparison with experimental results 
basis-set limits) is also made. 



27| (or computational 



A. Ground-State Energies and Ionization Energies at equilibrium geometry 

We first concentrate on calculating ground-state energies and ionization energies in equi- 
librium for a set of diatomic molecules with a singlet ground state. Calculations were per- 
formed for a number of different separation distances around the approximate equilibrium 
distance, after which a third-order polynomal was fitted to find the true energy minimum 
and equilibrium distance. The results calculated in a cc-pVDZ basis are presented in Table [B 

The ground-state energies for the molecules H2 to BeH2 show little difference (at most 4 
mH) between ADC(3) and FRPA. The differences for N2 and CO, containing double-triple 
bonds, are somewhat larger, of the order of 10 mH. The FRPAc ground-state energies tend 
to be close to the CCSD(T) results with a maximum deviation of 12 mH in case of BF. In 
general, ADC(3) deviates more from CCSD(T). 

The equilibrium bond distances show a larger spreading. The equilibrium bond distances 
for ADC(3) and FRPA have comparable deviations from the experimental values, and in the 
majority of cases are closer to the experimental value than the CCSD(T) results. The same 
conclusion can be made for the ionization energies, for which ADC(3) and FRPA outperform 
the coupled cluster results, when the experimental value is available. 

One remarkable fact is the lack of an equilibrium distance (no energy minimum) for N 2 
in both the ADC(3) and FRPA calculations without incorporating self-consistency at the 
level of the Hartree-Fock-like diagram. This example stresses the importance of a consistent 
treatment of the static self-energy. The inclusion of self-consistency in the calculations tends 
to adjust the results toward experiment, where needed. 

In order to compare with earlier ADC(3) calculations, we calculated vertical ionization 
energies for three diatomic molecules with the settings used in Ref. [28], i.e. at the ex- 
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perimental bond length and with the aug-cc-pVDZ basis set. The results are presented in 
Table HT1 T he present FTDAc results are in good agreement with the Dyson ADC(3) results 
in Ref . 28j . The differences are less than 2 mH and should probably be ascribed to a slightly 
different treatment of the HF-like self-energy. Compared to experiment, the mean absolute 
error is of the same order of magnitude for ADC(3) and FRPA. Note that there is a large 
deviation for the 2cr u -level of N2 in the FRPA which has a substantial influence on the mean 
error value. Apart from this level the mean absolute error of FTDAc and FRPAc is the 
same. 

We also checked the basis-set dependency of the results in Tables HHITl by performing cal- 
culations for HF in the cc-pVDZ, cc-pVTZ, aug-cc-pVDZ and aug-cc-pVTZ basis sets. The 
differences in ionization energies between DZ and TZ in Table HTT1 are of the order of 25 mH 
for the non-augmented and 10 mH for the augmented basis sets. The convergence behavior 
of the ground-state energies calculated with FRPAc are very comparable to CCSD(T). The 
weaker convergence in FRPA again demonstrates the importance of self-consistency for the 
Hartree-Fock-like diagram. 



B. Dissociation problems for H 2 

The FRPA fails to describe the correct dissociation behavior of diatomic molecules due 
to the appearance of instabilities in the RPA. The HF ground state becomes unstable with 
respect to ph-excitations in the dissociation limit. The RPA hamiltonian matrix is no longer 
positive-definite which results in complex solutions to the RPA equations. This is easily seen 
by analyzing H2 in a minimal basis set. The spatial wave functions are Is functions centered 
on the H- atoms A and B. These can be put in a bonding and anti-bonding combination, 
which will be the Hartree-Fock hole and particle state 

\b) = ^[\A) + \B)] 

\a) = j= 2 [\A)-\B)}. (39) 

These states are normalized to unity at great separation, which is the case we are interested 
in. 

The Hartree-Fock ground state is always the spatially symmetric state with positive parity 
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and spin: 



< F > = \bb)j,[\n)-\m (40) 

= «ULl° > • ( 41 ) 



The possible ph-excitations can only be formed by removing a bonding state and replacing 
it with an anti-bonding state. This results in a spin singlet and triplet. The energy of the 
triplet state 



1$ >= 



5=1 



al®a\ |0> (42) 

J M s 



is found to be 

< $|#|$ >= (ab\H\ab) - (ab\H\ba). (43) 

In the dissociation limit, the overlap and interaction matrix elements between the Is wave- 
functions for hydrogens A and B vanish, and the energy of the triplet state simply becomes 

< $\H\$ >= (AB\H\AB) + (BA\H\BA). (44) 

This is exactly the energy one would expect for the dissociation state where each hydrogen 
atom receives one electron, which is the exact ground state when the two hydrogens are 
separated by a large distance. The energy of this state is thus the exact ground-state 
energy, and automatically lower than the Hartree-Fock ground-state energy in this limit. As 
a result, a negative phonon energy occurs for this triplet state in the ph-TDA. In ph-RPA, 
the same mechanism gives rise to a complex phonon energy. 

This behavior is actually found both in the minimal basis set model for H 2 as in more 
realistic calculations. As an example of this behavior we have plotted the ground state 



energy for H 2 calculated in the cc-pVDZ basis set in Figure 4(a) For distances larger 



than approximately 1.2 A ph-RPA becomes unstable. At this distance the lowest ph-RPA 



eigenvalue in the spin-1 channel becomes zero, as can be seen in Figure 4(b) Beyond this 
distance ph-RPA acquires a complex eigenvalue. The ph-TDA eigenvalue becomes negative 
as well, which is unphysical for an excitation energy but does not pose any computational 
problems. A possible solution is to substitute the problematic spin-1 ph-RPA channel with 
its ph-TDA counterpart. The procedure then remains stable but is almost identical to 
FTDA. In any case, the use of TDA phonons does not guarantee a correct dissociation 
limit. Both the mixed and the pure TDA Faddeev method deviate substantially from the 
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FIG. 4. Demonstration of the problems in the dissociation limit for H2 in the cc-pVDZ basis set. 



(a) shows the ground-state energy for H2 calculated with the FRPA (solid line), FTDA (dashed 
line) and the mixed procedure (crosses) where the spin-1 channel of the RPA phonons has been 
replaced with TDA phonons compared to the exact full-CI result (dotted line). |(b)| shows the lowest 
ph-RPA and ph-TDA excitation energy in the spin-1 channel for H2 as a function of internuclear 
distance. 



exact full-configuration- interaction results in the dissociation limit. A more fundamental 
solution to this problem would probably be the fully self-consistent approach, where the 



propagator is allowed to have fragmented spectral strength 



This, however, implies a 



huge computational effort which lies beyond the scope of this paper. 



IV. CONCLUSION 

In this work we have investigated the application of the FRPA technique to small 
molecules. The computational cost of this method is not much higher than that of the more 
established ADC(3) method and in any case lower than the cost for CCSD(T). The results 
at equilibrium geometry are comparable in accuracy to the ones obtained with the ADC(3) 
method which is in line with the earlier atomic calculations. 

The self-consistent treatment of the Hartree-Fock diagram has a positive effect on the 
numerical results and should always be included. The possibility of complex eigenvalues in 
the RPA and FRPA eigenvalue equations is a problem that has to be kept in mind. We 
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have shown that RPA instabilities are bound to occur in the dissociation limit, when the 
Hartree-Fock propagator is used as a starting point. A possible way out is to increase the 
self-consistency by allowing propagators with fragmented single-particle strength, which will 
be the object of future research. 
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Appendix A: Derivation of the FRPA equations for i=3 

The product of the forward propagating uncorrelated 2plh propagator and the interaction 
vertex is needed to find an expression in function of RPA-amplitudes and the two-particle 
interaction. We will do this for the case % — 3, the other two cases are equivalent, but 
involving the T ph instead of the T pp . The combination of the free 2plh propagator and the 
vertex function can be written as 

[G (0)> r (s W W = \ f H / H E c'T («0 cgr (a - a) 

per 
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Here the explicit expression for the phonon propagator is needed. 
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The A pp are introduced as the product between the interaction and the normal RPA- 
amplitudes X pp and y pp . Due to the RPA-equations ( JT5]) . this correspondence can also 
be expressed as 

■ypp 

APP+ _ a P' m ( \o\ 



A 



V2 (eg* -e a - e fi ) 

3T 



pp— " Q/3,m 



V2 (eJT -e a - e p ) ' 

where the factor arises from the normalization condition for the pp RPA amplitudes and 
is not needed in case of ph RPA. 

After performing the necessary integrations over the intermediate energies, one arrives at 
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where in the second transition the property A (e Q + e^) = was used to simplify the 
relation. 



Appendix B: ADC(3) as special case of FRPA 

To show that ADC(3) is incorporated in FRPA one has to change the RPA-interactions 
with TDA-interactions. This can be done by setting the off-diagonal blocks in equations ( I15j) 
and (fT9|) to zero. As a result there are no backward propagating amplitudes y. The FRPA- 
equation fl33|) simplifies due to the disappearance of the HH^ . After projecting out the 

1(3 



spurious solutions we get the equation 

e Fd X = - (11 - 1U) ( U^D^U®- 1 + 2U®TW j (U - t ex ) X. (Bl) 

^ \i=l,2,3 / 

As an example we will again work out the term for i = 3 for the 2plh energies 

+25 7 , £ A*gJ (ef + - e Q - C/J ) A*£. (B2) 
By eliminating the TDA eigenvalues using their generating equations 

trxffn = ^ + e p ) {XZl ~ X&) + \ £ VafiMXfc (B3) 

A, /i 

and using the orthonormality of the TDA eigenvectors 

£ ^ PP+ = ^" A ^ - ^ a ^^ (B4) 

n 

we arrive at 

(^(3)^(3)^(3)-! + 2U^T^) a ^ = [(^ - $av$f3\) + 6, - 6,) + 3^, V ] . 

(B5) 

Similar steps have to be taken for the other two channels. The sum of the three channels 
after anti-symmetrization becomes 



- (11 - 1U) ( £ U^D^U {i) ~ l + 2rj«T« t ) (11 - ll e 

\i=l,2,3 / 

= <5 71 , (SaxSpfj, — S afM 5(sx) (e Q + e/j — e 7 ) 

This is exactly the same expression as in ADC (3). The Faddeev Tamm Dancoff Approxi- 
mation (FTDA) and ADC(3) are completely equivalent. 
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TABLE I. FRPA results for some diatomic molecules and BeH2 in a cc-pVDZ basis set. The 
ground-state energy Eq and vertical ionization energy I are in Hartree, equilibrium bond distance 
ro is in Angstrom. FRPA and FTDA refer to the calculations after the first iteration, while FRPAc 
and FTDAc refer to the calculations where consistency on the Hartree-Fock level was applied. The 
calculated data are compared to the high-level ab-initio method CCSD(T) where available and to 

n 

experimental data or exact calculations from Ref. [271 ] . 



Molecule 




FTDA 


FTDAc 


FRPA 


FRPAc 


CCSD(T) 


Expt. 


H 2 


E 


-1.170 


-1.161 


-1.170 


-1.161 


-1.164 


-1.175 




ro 


0.769 


0.757 


0.770 


0.757 


0.761 


0.741 




I 


0.594 


0.589 


0.594 


0.589 


0.583 


0.591 


HF 


Eo 


-100.175 


-100.224 


-100.173 


-100.228 


-100.228 


- 




ro 


0.904 


0.916 


0.897 


0.913 


0.920 


0.917 




I 


0.577 


0.577 


0.572 


0.571 


0.628 


0.592 


HC1 


Eo 


-460.295 


-460.256 


-460.293 


-460.258 


-460.254 


- 




ro 


1.314 


1.297 


1.314 


1.293 


1.290 


1.275 




I 


0.457 


0.450 


0.457 


0.450 


0.471 


- 


BF 


E 


-124.331 


-124.365 


-124.332 


-124.368 


-124.380 


- 




ro 


1.285 


1.284 


1.305 


1.285 


1.295 


1.267 




I 


0.417 


0.395 


0.431 


0.402 


0.406 




BeH 2 


Eo 


-15.855 


-15.831 


-15.856 


-15.832 


-15.835 






r 


2.747 


2.674 


2.766 


2.674 


2.678 


2.680 




I 


0.437 


0.433 


0.435 


0.432 


0.446.1 




N 2 


Eo 




-109.258 




-109.272 


-109.276 






ro 




1.104 




1.106 


1.119 


1.098 




I 




0.565 




0.544 


0.602f. 


0.573 


CO 


Eo 


-113.096 


-113.037 


-113.100 


-113.048 


-113.055 






ro 


1.140 


1.130 


1.133 


1.123 


1.145 


1.128 




I 


0.529 


0.503 


0.523 


0.494 


0.550 a 


0.515 



Only up to CCD level 



20 



TABLE II. Vertical ionization energies in Hartree calculated in the aug-cc-pVDZ basis set. The 
values between braces are calculated without the la u -level of N2. Experimental values are from 
Ref. 



cc-pVDZ 



FRPA 

aug-cc-pVDZ cc-pVTZ 



Expt. 



aug-cc-pVTZ 



Molecule 



Level 



HF 



l7T 

3a 



15.46 
19.57 



16.06 
20.01 



16.18 
20.06 



16.33 
20.21 



16.05 
20.00 



FRPAc Expt. 

cc-pVDZ aug-cc-pVDZ cc-pVTZ aug-cc-pVTZ 

Molecule Level 
HF 

lvr 15.53 16.34 16.17 16.42 16.05 

3a 19.54 20.24 20.00 20.27 20.00 
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TABLE III. Ground state energies and vertical ionization energies in Hartree for HF, calculated in 

f. In- 



different basis sets. Experimental values are from Ref. 



28|, CCSD(T) values are from Ref. 



Molecule 


Level 


HF 


FTDA 


FTDAc 


FRPA 


FRPAc 


Expt. 


HF 


















Itt 


0.651 


0.596 


0.605 


0.590 


0.601 


0.592 




3a 


0.771 


0.740 


0.747 


0.736 


0.744 


0.735 


CO 


















5a 


0.555 


0.532 


0.510 


0.528 


0.503 


0.515 




Itt 


0.641 


0.626 


0.622 


0.623 


0.619 


0.621 




4a 


0.808 


0.737 


0.739 


0.715 


0.720 


0.724 


N 2 


















3 a „ 


0.634 


0.593 


0.575 


0.579 


0.558 


0.573 






0.615 


0.632 


0.618 


0.651 


0.630 


0.624 




2a u 


0.781 


0.711 


0.698 


0.672 


0.658 


0.690 




A abs (mH) 


49 (44) 


12 (10) 


8(8) 


10 (9) 


11 (8) 






A max (mH) 


91 (84) 


21 (20) 


15 (15) 


27 (27) 


32 (15) 
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